Mixed-spin cluster expansion for a quasi-one-dimensional Haldane system 
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We present a novel mixed-spin cluster expansion method 
for a quasi-one-dimensional Haldane system with bond al- 
ternation. By mapping the s — 1 antiferromagnetic spin 
model on square and cubic lattices to the equivalent mixed- 
spin model, we study the competition among the Haldane, the 
dimer and the magnetically ordered phases. The mixed-spin 
cluster expansion proposed here allows us to directly deal with 
the Haldane phase, which may not be reached by standard 
series expansion methods. The phase diagram is determined 
rather precisely by making use of an additional symmetry 
property in the effective mixed-spin model introduced. 

Low-dimensional spin systems with the spin gap for 
the excitation spectrum have-been extensively studied 
since the Haldane conjecture^! which clarified that the 
gap formation in the integer-spin Heisenberg chain re- 
flects the topological nature of spins. Recent extensive 
experimental and theoretical investigations on the stabil- 
ity of the Haldane system against various perturbations 
have been providing a variety of interesting topics. The 
instability of the spin-gap phase in the s = 1 spin mod- 
els has been studied in detail so far for one-dimensional 
(ID) systems. For instance, the effect of the bond alter- 
nation is understood qualitatively well by the non-linear 
sigma modela, as well as the valence bond solid (VBS) 
approach.!] The accurate critical point between the dimer 
and the Haldanejphases has been further obtained by the 
series expansion,!!! the exact diagonalizationJ§El the quan- 
tum Monte Carlo simulations!] rand the density matrix 
renormalization group (DMRG) J3 On the other hand, the 
s = 1 spin systems with the 2D or 3D structures have not 
been studied so well, although the effects of the antifer- 
romagnetic correlations due to the interchain couplings 
should be important for real materials. So far, Sakai and 
Takahashila investigated a quasi-lD s = 1 spin system by 
combining the mean field theory with the exact diago- 
nalization results for the spin chain, and gave a rough 
estimate for the phase-transition point to the antiferro- 
magnetic phase. 

In this paper, we systematically study how the Hal- 
dane and the dimer phases for the s = 1 antiferromag- 
netic chain are driven to the magnetically ordered phase 
in 2D and 3D systems by exploiting the series expan- 
sion techniques. In particular, we propose a mixed-spin 
cluster expansion by mapping the s — 1 spin model to 
the equivalent mixed-spin model, which allows us to deal 
with the Haldane phase. This new approach is a realiza- 
tion of the notion of the VBS in a perturbation theory. 
We determine the phase diagram rather precisely both 
for the 2D and 3D cases by computing the spin excita- 
tion gap and the staggered susceptibility. 



Let us first consider the s = 1 antiferromagnetic quan- 
tum spin system on a 2D square lattice, which is de- 
scribed by the Hamiltonian 



JS 
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(1) 



where J is the interchain coupling and Sjj is the s = 1 
operator at the (i,j)-th site in the (x — y) plane. Here we 
have introduced the bond-alternation parameter a(0 < 
a < 1) along the x direction, Tj = 1(a) for even (odd) i. 
All the exchange couplings are assumed to be antiferro- 
magnetic. 

We employ the series expansion method developed by 
Singh, Gelfand and HuseB Since this method combines 
the conventional perturbation theory with the cluster ex- 
pansion, it has an advantage to deal with the spin system 
in higher dimensions even for the cases for which the reli- 
able results are difficult to be obtained by the exact diag- 
onalization, the DMRG, etc. In fact, the series expansion 
method has been successfully-applied to the 2DpSpin sys- 
tems with. various structures, ErE9 Kondo lattice,til bilayer 
systems,L3 etc. However, to apply the series expansion 
technique to the present system including the Haldane 
phase, a nontrivial generalization is needed, since a naive 
cluster expansion may not describe the Haldane state. 
For instance, the dimer state is adiabatically connected 
to the isolated s = 1 dimers, but the Haldane state does 
not have its analogue in the isolated local singlets com- 
posed of several s = 1 spins. To overcome this problem, 
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FIG. 1. (a) The s = 1 spin chain, which is decomposed 
into (b) the mixed-spin chain in the Haldane phase and into 
(c) that in the dimer phase. Large and small solid circles 
represent the s = 1 spin and s = 1/2 spins, respectively. The 
solid lines in (b) and (c) indicate the strong bonds which make 
local singlets whereas the dashed lines the weak bonds which 
are treated perturbatively. 

we wish to recall the notion of the VBsi, which captures 
the essence of the Haldane-gap formation. To realize this 
idea in the series expansion, we first divide half of the 
s = 1 spins into two s = 1/2 spins as schematically shown 
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in Fig. pi3 and map the system to the mixed-spin sys- 
tem which is equivalent to the original model except for 
a trivial isolated excited mode. As a starting configura- 
tion in the perturbative expansion, we can then consider 
two types of the mixed-spin cluster singlets formed by 
the solid lines in Figs. |l| (b) and (c). It is seen that by 
starting from the configuration (b) we can directly deal 
with the Haldane phase since it has the structure of the 
Haldane state in the VBS picture, whereas if the config- 
uration (c) is chosen, we naturally end up with the stan- 
dard dimer expansion. The above mapping thus gives us 
an important message that the Haldane phase is adiabat- 
ically connected to the isolated mixed-spin singlet states 
in Fig. \^(b), and thereby can be treated by the mixed- 
spin cluster expansion method. The resulting cluster ex- 
pansion around the isolated mixed-spin singlets should 
provide a quite powerful method, which enables us to 
deal with the competition among the Haldane phase, the 
dimer phase and the magnetically ordered phase in 2D 
and 3D systems. 

Let us begin with the quantum phase transition be- 
tween the Haldane phase and the antiferromagnetic 
phase in 2D s — 1 spin system with bond alternation. 
To this end, we consider the effective 2D mixed-spin sys- 
tem shown in Fig. |[ In this figure, the large (small) 




FIG. 2. (a) The 2D s — 1 spin model and (b) the corre- 
sponding mixed-spin system. 

circle represents the s = l(s = 1/2) spin. The bold 
solid, the thin solid and the dashed lines indicate the 
coupling constant 1, A and JA, respectively. In this fig- 
ure, the model without bond alternation is drawn for 
simplicity. We note that the mixed-spin system repro- 
duces the original 2D spin system at A = 1. To per- 
form the cluster expansion, the Hamiltonian is divided 
into two parts as H = ^2 Sj • Sj r + A 53 S< • Sj . The first 
term is the unperturbed Hamiltonian which stabilizes the 
isolated mixed-spin cluster singlets. The corresponding 
mixed-spin cluster has the configuration, 1/2 o 1 o 1/2, 
which is formed by the antiferromagnetic couplings 1 and 



a. These isolated clusters have the singlet ground state 
with the spin gap A = (3a + 3 - y/9 - 14a + 9a 2 )/4. 
The perturbed part of the Hamiltonian labeled by A con- 
nects these isolated mixed-spin singlets to form a 2D net- 
work and thus enhances the antiferromagnetic correla- 
tion. We compute the staggered susceptibility xaf, and 
the singlet-triplet excitation gap A at the ordering wave 
vector. These quantities are then expanded as a power 
series in A. We finally determine the phase boundary 
by the divergent staggered susceptibility and the vanish- 
ing spin gap, Atthich are estimated by applying the Pade 
approximantai-3 to the quantities obtained up to the finite 
order in A. 

To confirm how well our mixed-spin cluster approach 
works, we first investigate the s = 1 spin chain without 
bond alternation. Performing the mixed-spin cluster ex- 
pansion, we calculate the ground state energy E g , the 
staggered susceptibility xaf and the singlet-triplet exci- 
tation gap A up to the eleventh, the fifth and the sev- 
enth order, respectively. At first sight, the order in the 
series for the staggered susceptibility and the excitation 
gap might not be high enough to produce the accurate 
values at A = 1 (the Haldane_.point) by means of the 
ordinary differential methods.EJ It is remarkable, how- 
ever, that there exists an additional symmetry property 
like Q(\) = XQ(1/X) for each quantity Q in our effec- 
tive mixed-spin chain, which enables us to expand the 
quantity Q as a power series even around A = 1. Fit- 
ting this power series with that obtained by the cluster 
expansion, we end up with the rather accurate values, 
Eg = -1.4022, xaf = 19.6 and A = 0.404, which are 
compared with those of the Monte Carla simulations: 
E g = -1.4015 ± 0.0005, A = 0.41 in ref.,E>l and also the 
exact dijagonalization A = 0.411±0.001, xaf = 18.4±1.3 
in ref.E3 
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FIG. 3. Phase diagram for the 2D s = 1 quantum spin sys- 
tem with bond alternation a. The phase boundary between 
the Haldane phase and the ordered phase is determined by the 
mixed-spin cluster expansion. The left solid (the left dashed) 
phase boundary around the dimer phase is determined by the 
dimer expansion with the biased [2/3] Pade approximants for 
the excitation gap (the staggered susceptibility). 

To study the quantum phase transition on a 2D lattice 
by increasing the interchain couplings, we evaluate the 
singlet-triplet excitation gap A by means of the mixed- 
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spin cluster expansion up to the fifth order in A for vari- 
ous choices of a and J. It is sufficient to consider the pa- 
rameter regime near A = 1 to discuss the original Haldane 
system. In the case without bond alternation (a = 1), 
applying the Dlog Pade approximants to the spin gap, the 
critical value J c = 0.056 ± 0.001 and the critical expo- 
nent v = 1.86 ± 0.08 are obtained. Our results for a = 1 
are much more accurate than those of the mean field 
theory combined with the exact diagonalizationja which 
claimed the critical value to be J c > 0.025. We here 
note that the obtained critical exponent is different from 
the value v p_0.7lll3 expected for the 3D classical Heisen- 
berg model.ES This implies that the quantum phase tran- 
sition in our generalized mixed-spin model does not be- 
long to the universality class of the 3D classical Heisen- 
berg model in generic cases, although it should do in the 
specific case A = 1. Assuming that the spin gap in the 
vicinity of the transition point vanishes with the same 
exponent even for the Haldane system with bond alter- 
nation, we determine the phase boundary shown as the 
dots with error bars in Fig. ^. The error bars come from 
the different values obtained by different biased Pade ap- 
proximants employed: [1/2], [2/1], [2/2], [2/3], [3/2] ap- 
proximants. Since the error bars increase with the de- 
crease of a away from unity, it seems difficult to deter- 
mine the phase boundary in the region close to the dimer 
phase. However, it is to be noted that this phase diagram 
should have the symmetry property as J {a) — aJ(l/ot). 
Taking this into account, we can thus determine rather 
precisely the phase boundary between the Haldane phase 
and the antiferromagnetic phase, which is drawn by the 
solid line in Fig. ||[ We shall see momentarily that the 
critical point between the dimer and the Haldane phases 
determined in this procedure is quite consistent with that 
obtained by the dimer expansion. 

Let us now turn to the dimer phase. In this case, our 
mixed-spin cluster expansion is equivalent to the stan- 
dard dimer expansionEa We perform the dimer expan- 
sion of the staggered susceptibility and the spin gap up 
to the fifth and the sixth order in A for various J, respec- 
tively. To estimate the phase boundary which separates 
the dimer phase and the antiferromagnetic phase, we use 
the ordinary Pade approximantsEj as well as the biased 
Pade approximants, for which the phase transition is as- 
sumed to belong to thfi-.unversality class of the 3D classi- 
cal Heisenberg model. La Using these Pade approximants, 
we arrive at the phase diagram shown in Fig. ||. When 
J = with small a, the system is reduced to the iso- 
lated s = 1 bond-alternating chain, which is known to 
have disordered ground state with the spin gap due to 
the dimer singlet. Increasing the parameter J and a, the 
antiferromagnetic correlation grows up, and the quantum 
phase transition to the magnetically ordered state occurs. 
We wish to note that the critical point (a, J) = (0.59, 0), 
which is determined from the series expansion of the spin 
gap, separates the Haldane phase, the dimer phase and 
the antiferromagnetically ordered phase in Fig. |3|. Since 
the system in this case is reduced to the independent 



s = 1 spin chains with bond alternation, our numerical 
results reproduce the well-known factQLl that the ground 
state of the reduced chain with a c — 0.59 is in a critical 
phase with neither the spin gap nor the long-range order. 
To confirm how accurate our results for 2D cases are, we 
have directly analyzed the spin chain (J = 0) by applying 
the Dlog Pade approximants to the spin gap computed 
up to the eighth order. This gives a c = 0.612 ± 0.004, 
which is close to the value 0.59 obtained above, and also 
to 0.60 ± 0.01 obtained by DMRG.I Judging from these 
results, we can say that our phase boundary determined 
by the excitation gap in Fig. [| is quite accurate, while 
that by the staggered susceptibility has a slight deviation 
only around the critical point. 
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FIG. 4. Dispersion relations of the spin-triplet excited 
states for the s = 1 chain with bond alternation a. The 
results for a = 0,0.3 and a c (= 0.59) are obtained by the 
dimer expansion, while the other for the Haldane phase is ob- 
tained by the mixed-spin cluster expansion with the help of a 
symmetry property (see the text). 

In order to demonstrate that our approach is also 
powerful to compute the elementary excitation with fi- 
nite momenta, we show the calculated dispersion rela- 
tion in Fig. |^ along the specific line of J = 0. Reflect- 
ing the isolated spin-singlet structure, the Brillouin zone 
becomes half of the original one. In the dimer phase 
< a < a c , using the first order inhomogeneous differ- 
ential method,E-3 we can obtain the dispersion relation. 
Here, to obtain the dispersion for the Haldane phase, we 
have again made use of the additional symmetry prop- 
erty inherent in the effective mixed-spin model mentioned 
above. It is to be emphasized that such a precise dis- 
persion is obtained within the lower-order perturbations, 
which is indeed due to the additional symmetry we have 
used. 

We now move to the 3D system. The advantage of our 
approach is particularly remarkable for the 3D problem 
because other numerical methods may often meet some 
difficulties to treat a large system in the 3D case. We here 
consider a cubic lattice system by adding the interchain 
couplings J in the z direction to the spin model discussed 
above. By extending our treatment to the 3D system, 
we thus study the competition between the two kind of 
gapped states and the antiferromagnetic state. Applying 
the dimer expansion to calculate the spin gap and the 
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staggered susceptibility up to the fifth order and using 
the Dlog [2/2] Pade approximants, we first determine the 
phase boundary which separates the dimer and the anti- 
ferromagnetic ordered phase in Fig. ||. When a = 0, our 
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FIG. 5. Phase diagram for the 3D s = 1 spin system with 
bond alternation a. 

system reproduces the s = 1 bilayer Heisenberg model. 
Increasing the inter-dimer coupling J from zero, the an- 
tiferromagnetic correlation grows up and the quantum 
phase transition occurs at J c = 0.143 ± 0.006. We note 
that the quantum phase transitions in the biiwer model 
have already been studied by Gelfand et alE2l with the 
series expansion method. On the other hand, to observe 
the phase transition from the Haldane phase to the or- 
dered phase, we further perform the mixed-spin cluster 
expansion up to the fourth order for both of the above 
quantities. In the homogeneous case (a = 1), by analyz- 
ing the data in terms of various Dlog Pade approximants 
we end up with the critical point J c = 0.026 ± 0.001, 
which is consistent with those of the non-linear a model 
approachtll and the. mean field theory combined with the 
numerical method.El The phase diagram thus determined 
is shown in Fig. |^. 

In summary, we have investigated the quantum phase 
transitions for the s — 1 quantum systems with the 2D 
and 3D structures. Using the series expansion, we have 
discussed how the dimer phase and the Haldane phase 
realized in ID compete with the magnetically ordered 
phase in higher dimensions. In particular, we have pro- 
posed a novel approach based on the mixed-spin cluster 
expansion which realizes the idea of the VBS in the per- 
turbation theory. This new approach has made it pos- 
sible to treat the Haldane phase in the series expansion 
framework, which was not dealt with so far by ordinary 
series expansion methods. For the spin chain case, we 
have obtained faipkj-good results comparable to other nu- 
merical methods.E£rE3 For the 2D and 3D cases, the phase 
diagram has been determined rather precisely by making 
use of an additional symmetry property in the effective 
mixed-spin model. It is quite interesting to futher apply 
the mixed-spin cluster approach to the frustrated case, 
the anisotropic case, etc., in quasi-lD Haldane systems, 
which is now under consideration. 
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